Mean-field approximation for a limit order driven market model 
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investigated |L£ 

The above mentioned models mostly do not take 
into account the realistic details of the price formation 
through the book of orders. This mechanism was imple- 
mented in the model set up by Maslov |35| and a similar 
perspective was then used in a recent series of papers by 
Matassini and Franci 136-38 . The book of orders per- 
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The mean-field variant of the model of limit order driven market introduced recently by Maslov 
is formulated and solved. The agents do not have any strategies and the memory of the system 
is kept within the order book. We show that he evolution of the order book is governed by a 
matrix multiplicative process. The resulting stationary distribution of step-to-step price changes is 
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I. INTRODUCTION 

The complexity of market behavior, seen as a particu- 
lar example of a natural phenomenon, fascinated physi- 
cists since long time ago . The main source of interest 
comes from a kind of critical behavior, made explicit by 
the power-law distribution and scaling in the economic 
time series, first observed by Mandelbrot (see [|| and ref- 
erences therein) and studied in detail by Mantegna and 
Stanley |3|-|5| and subsequently by many others (see e. g. 
10)- 

The scaling and multifractal properties call for an ex- 
planation in terms of a model mimicking the behavior 
of individual agents in the market. The physicist's opti- 
mism in looking for such a model might be strengthened 
by recent successes in modeling other social phenomena, 
e. g. in the cellular-automata models of traffic The 
idea consists in assuming that all the overwhelming com- 
plexity of a human being is irrelevant in certain special 
conditions: when driving a car, only a very basic set 
of behaviors is at work. Similarly, it is assumed that a 
trading agent, when put on the floor, follows only a lim- 
ited set of instincts or acquired patterns. Therefore, in 
this approach the economic complexity is not due to the 
intrinsic complexity of each agent (as usual hand-waving 
arguments by liberal opponents of "reductionism" state), 
but an emergent property of a large set of non-linearly 
interacting simple units. 

Many microscopic stock market models emerged dur- 
ing last several years. One of the first ones was the 
model introduced by Levy, Levy and Solomon 13 1^] 
which captures essential features of the price fluctuations 
and explains also the power law distribution of investor's 
wealth, which is the famous Pareto law. 

Another approach was used in the model of Bak, 
Paczuski and Shubik Jl6| |: buyers and sellers are repre- 
sented by particles subject to a reaction-diffusion process. 
Introduction of a non-trivial strategy of the agents lead 
to a realistic value of the Hurst exponent for the price 
fluctuations. The simpler version of the model was then 
solved analytically |17| . Variety of other approaches were 



spective to market modeling was empirically investigated 
in §9)0. 

With such a diversity of models, most of which giv- 
ing plausible explanation of observed facts a question 
arises, whether there is a common mechanism behind 
various approaches, making them essentially equiva- 
lent. Indeed, it was found that such a mechanism may 
be the multiplicative stochastic process repelled from 
zero, or multiplicative-additive process. It was studied 
thoroughly by various authors and in diverse contexts 
fll5] , [n| -|5l| . The goal we pose in this work is to show 
that essentially the same mechanism is responsible for 
the power-law distribution of price changes also in the 
limit-order model. 



II. LIMIT ORDER DRIVEN MARKET MODEL 

Recently Maslov |35|,[59| proposed a model, based on 
the assumption that there are two kinds of market partic- 
ipants. Prudent investors place their orders at prescribed 
price and a trade occurs as soon as there is anybody ac- 
cepting that price. On the other hand, speculators buy 
and sell at any moment at the price which is available 
in the market. The price signal p(t) was found to have 
power-law spectrum, with Hurst exponent H = 1/4. The 
price changes during unit time interval x — p(t+ 1) — p(t) 
have probability distribution which follow clear power 
law P{x) ~ .t _( - 1+q ) in two regimes. For small x the ex- 
ponent is 1 + ol\ = 0.6 ±0.1, while in the regime of large 
price changes the exponent is 1 + «2 = 3 ± 0.2. These 
values are to be compared to the experimentally found 
1 + ax ~ 2.5 and 1 + a 2 ~ 4 respectively 

The model [B5f can be described as follows. There are 



1 



orders to buy and sell placed on a straight line, which is 
the axis of the price x. In stable situation, all buy orders 
are lower than all sell orders, so that we can describe the 
state by single function p(x), density of the orders, and a 
number £ which is the last realized price. Then, all x < £ 
correspond to buy, all x > £ to sell orders. 

Two events can change the state. First, new limit or- 
ders may be dropped, such that p(x) — > p(x) + r\{x — £) 
We suppose the function rj(x) equal in all events and sym- 
metric, r](x) = 7](—x). Second, a market-price order can 
arrive. An order to buy an amount s results to clearing 
all sell orders up to the price £ + x + , where 

J p(x)dx = s (1) 

The new price is then £ — > £ + so that x + is the price 
increment, while the new density is p(x) —> (1 — 9(x — 
+ x + — x j) p(x). Analogical formulae will hold for 
the sell order. 

As we can see, there are no strategies which would 
lead the agents to perform specific actions. The model is 
barely stochastic. The long-term memory of the system 
and thus a possible power-law behavior stems from the 
order book, or the time-dependent density function p(x) 
which may keep arbitrarily old orders. 

III. MATRIX FORMULATION 

Our essential approximation to this model will consist 
in supposing uniform density of orders on each side from 
the current price level £. In reality, both dropping new 
limit orders and clearing them by market orders makes 
the density of states uneven and fluctuating. When sup- 
posing that after an event the uniform density of states is 
restored, we make a kind of "mean-field" approximation: 
the actual position of each limit order is not important, 
as if they were freely moving particles making an effective 
medium, within which the price fluctuates. High density 
of the medium will result in smaller price fluctuations 
and vice versa. 

The density on the upper side will be denoted p + , on 
the lower side p_ . It is convenient to describe the densi- 
ties in terms of the potential price changes, which would 
occur if a market-price order arrives. They are simply 
x + = s/p + for buy and x_ = s/p- for sell order. The 

numbers x± form a vector X = ^ ^ + ^ which performs 

a stochastic process, as the densities p± and therefore 
the numbers x± are updated after arrival of each order. 
The dynamical rules of the process represent a simplified 
version of the limit order driven dynamics. 

There are be three types of events: (i) dropping of limit 
orders (ii) market-price order to buy, and (iii) market- 
price order to sell. We suppose that all market price 
orders have the same volume s and all limit-order events 



the same volume v. Further we assume that market or- 
ders to buy and sell occur with the same probability. In 
order to keep the total number of limit orders constant on 
average, we should suppose that at given moment there 
is probability p = s/(s + v) to drop a limit orders. Each 
of the market-price events (ii) and (iii) have then equal 
probability (I — p)/2. (Here we tacitly assume that the 
limit orders which are not met last forever. One can also 
investigate more realistic variants, where the limit orders 
slowly die out). 

Let us investigate first the consequence of an arrival 
of a market-price order to buy. By definition, the price 
level increases by x+ . As the density of orders is constant, 
average density on the right-hand side from the new price 
is unchanged. If now another buy order arrives, it finds 
the same density and the price change is the same too. 
Therefore, the new value of x + is equal to the old one, 
x + — > x + . On the other hand, if now a new market order 
to sell arrives, there are no limit orders in the interval 
of width x+ below the current price level, and when we 
go further down, there is a constant density s/x-. As a 
result, the price decreases by x + + X-. Hence, the new 
value of is — > x + + X-. 

To sum it up, the effect of the buy order consists in 
the replacement 

x+ ► ,„\ 

X- — > X + + X- . ^ ' 

It can be expressed in matrix form 

X -> X' = T+X (3) 

where 




Similarly, for the action of a sell order we get 

X -> X' = T-X (5) 

where 

r-(il). 

Now we turn to the changes due to dropping limit or- 
ders. It is necessary to specify the function f}{x), rep- 
resenting the average volume of orders set at distance 
x from the current price. As we already mentioned, we 
suppose it to be an even function. Moreover, the vol- 
ume was supposed to be fixed, J r](x)dx = v. We apply 
here the simplest choice r}{x) = |(<5(x — d) + 5(x + d)), 
which means that all new orders are placed at the same 
distance d from the current price, either below (buy) or 
above (sell). This distribution reflects the fact, that the 
limit orders are not typically set arbitrarily close to the 
current price, but there is a certain minimum offset d. 
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Dropping limit orders affects the vector X according 
to the formula 



x ± -> -(3 - l/p)x± 



(7) 



Indeed, a buy (or sell) order will annihilate the amount 
v/2 from the just deposited limit order and amount s — 
v/2 from the original density of old limit orders. The 
shift is therefore x± = (s — v/2)/p±. Writing v in terms 
of the probability p, i. e. v = (1/p — l)s, we obtain the 
formula (0). So, in matrix form we have 



X ^ X' = SX 



where 



i o 

1 



(8) 



(9) 



The price changes only after a market order is is- 
sued, while dropping limit orders leaves the price un- 
changed. So, between two subsequent shifts of the price, 
to > limit orders can arrive, with probability P m {m) = 
(1 —p)p m . The change of the vector X due to one market 
order and to limit orders is X — > S m T± X. When calcu- 
lating the evolution of the probability distribution for X, 
we should sum over all possible realizations. Hence, the 
probability distribution for the vector X should satisfy 
the equation 

Px{X) = (10) 
2 E E / dX'P m (m)P x (X')5(X - S m T a X') 

<7=± m=0 J 

in the stationary state. 

IV. DISTRIBUTION OF PRICE CHANGES 

We will make a further approximation at this stage. 
The matrix S is simply a unit matrix multiplied by a 
constant. If the same were true also for the matrices T±, 
the process would be reduced to a simple multiplicative 
random walk, whose properties are well known and their 
relevance in modeling price fluctuations is testified by a 
series of models, as mentioned in the Introduction. 

Our approximation will consist first in replacing the 
matrices T± by the average T = \(T+ + T_) and fur- 
thermore, we will take only the highest eigenvalue of the 
matrix T, which is 3/2. Then, instead of a pair of price 
changes x+ and X- we have a single scalar quantity x, 
describing the absolute value of the price change. 

Note that the same results can be obtained by assum- 
ing from beginning, that x + = x_, i.e. that the density 
of states is equal on both sides of the price level. This 
means, that we make a further "mean-field" approxima- 
tion, suppressing not only the fluctuations along the price 



axis, but also fluctuations from one side to the other of 
the already averaged density of states. 

This way we define our multiplicative random process. 
The fact that there is a small but finite offset d in placing 
the limit orders ensures that the values of x± (therefore 
also of x) cannot be smaller than d. This feature plays the 
role of "repulsion from zero" , which was found essential 
for establishing the power- law tai ls BEsf and is usually 
guaranteed by the additive term p3^6]. 

For the probability distribution of the price changes we 
obtain 



(11) 




and assuming a power-law tail of the probability distri- 
bution in the form P(x) ~ x^ 1 ^ we obtain the following 
equation for the exponent 



m=0 



(1 -p)p 




(12) 



As can be easily checked, apart from the trivial solution 
a = it has a non-trivial solution a = 1, independent 
of p. Therefore, the distribution of price changes has a 
power-law tail 

-2 



P{x) 



(13) 



Note that the calculation could be further simplified by 
writing the equation analogical to ([ll]) , relating the prob- 
ability distribution just after single step. Then, instead 
of the sequence of steps consisting of one market order 
followed by to limit orders we have one step being ei- 
ther limit or market order. This equation gives precisely 
the same power-law tail. However, such an approach is 
slightly inconsistent, because setting a limit order does 
not imply any trade, thus the price change at this mo- 
ment is zero. 



V. PRICE AUTOCORRELATION FUNCTION 

One of the well-known facts about financial data se- 
ries is the negative short-time autocorrelation of price 
changes 0. Here we will show, how this effect naturally 
emerges from the matrix nature of our stochastic process. 

We will compute the autocorrelation function defined 

as 

(X{t)x(t + T)} 



C(t,t- 



(14) 



where x(t) is the actual price change at time t and r > 1. 
We will denote M(t) e {5,T + ,T_} the matrix describ- 
ing the action performed at time t and Pm(M) its prob- 
ability distribution. Of course, we introduced already 
P M (T_ ) = P M (T+) = (1 - p)/2 and P M (5) = p. 
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Now we introduce the function of taking the price 
change from the vector X — [ x+ ) 

X(X; M) = x+ if M = T+ 

= if M = S (15) 
= a;_ if M = T- . 

Note that the operator X is linear in the argument X. 
Then 

(x(t)x(t +t)> = JaXP x {X)C{X) (16) 

where 

T 

C(X)=J2- E I[PM(M(t + i))x (17) 

M(t) M(t+r) i=0 

x X{X\ M(t)) X(M(t + t- l)...M(t)X] M(t + r)) 

We find easily, using the linearity if the operator X 
and introducing the sign vector E = (1,-1) 

C{X) = 1^e(Y j P m {M)M^ x (18) 

x(x+T + - x_T_)A . 

The multiplication by E extracts only the lower one of 
the two eigenvalues of the averaged matrix 

5>.(M)M=I(j+^) . (19) 

The lower eigenvalue is p, hence 

C{X) = -{l-p)p r - 1 x + x- . (20) 

We can calculate similarly the corresponding expres- 
sion for the denominator of the equation ( |l5| ) . We obtain 
at the end 

^ ,+T >-^$F&- (21 » 

We can clearly observe the negative autocorrelation 
which decays with characteristic time which depends on 
the relative frequency of putting the limit and market 
orders, measured by the probability p. The result ob- 
tained suffer from the divergence of second moments of 
the variables x+, x_, resulting from the power- law tail 
calculated in the last section. However, only the ratio 
of the moments enter the formula (|2l]). Moreover, if we 
suppose as an initial condition a distribution for x+, X- 
with finite moments, the moments will remain finite for 
any finite time and we naturally expect that the ratio of 
the second moments will converge to a finite value even 
if the second moments themselves diverge. 



VI. CONCLUSIONS 

In conclusion, we solved in the mean-field approxima- 
tion the Maslov model of stock market fluctuations. We 
found a stationary distribution of price changes with a 
power law tail with the exponent 1 + a = 2, which is 
within the Levy stable region. We found negative short- 
time autocorrelation of the price changes, decaying ex- 
ponentially with time. The relaxation time depends on 
the relative frequency of putting market orders and limit 
orders: the decay is slower if a market order comes only 
after more limit orders. This is intuitively clear, because 
it is the market order that ensures the liquidity. 

Our result differs in two important points from the 
simulations of Maslov |55). First, the numerical value 
of the power-law tail exponent is different. This can be 
attributed to the approximation we made in the form of 
the density of orders p(x). Indeed, we assumed, as a zero 
approximation, constant density. On the other hand, it is 
known from the solution of the reaction-diffusion model 
of market JlTj that the density may have complicated 
non-trivial form. Another source of the difference may 
be the neglect of fluctuations. 

Another difference consist in lacking the second power- 
law regime for small price changes. However, as discussed 
in J35|, this second and different power law comes from 
the fact, that the new limit orders may be placed far- 
ther than the reach of the price change. As we implicitly 
supposed that the new orders are put to very small dis- 
tance d from the current price, we can observe only the 
distribution for price changes larger than d. 
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